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Abstract 

We consider the sensitivity, with respect to a parameter 9, of paramet- 
ric families of operators Ae, vectors ne corresponding to the adjoints Al 
of Ae via Agire = and one parameter semigroups t i-)- e*'*" . We display 
formulas relating weak differentiability of S i— > vre (at 6 = 0) to weak dif- 
ferentiability of 6 ^ AgTvo and [e'^''*]*7ro. We give two applications: The 
first one concerns the sensitivity of the Ornstein-Uhlenbeck process with 
respect to its location parameter. The second one provides new insights 
regarding the Wright-Fisher diffusion for small mutation parameter. 
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Introduction 



Sensitivities of parametric families of dynamical systems (with respect to the 
parameter) have been studied in the context of stochastic processes [TU] as well 
as partial differential equations [2] and are useful tools in optimization and con- 
trol. We consider sensitivities in the setting of one parameter semigroups. This 
setting constitutes a unifying approach to continuous time Markov processes 
[T] and linear PDEs [3], [T]. In particular it allows an elegant treatment of the 
sensitivities of the Wright -Fisher diffusion. 
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We consider one parameter semigroup^ U := (C^(i))te[o,oo), that can be 
represented as U{t) — e*"^ for some linear operator A : E ^ E on some linear 
space E. This does not exclude generators A that are usually considered as 
'unbounded', since we do in general not suppose that i? is a Banach space. (For 
information on the Hille-Yosida Generation Theorem providing U {t) = e*^ for 
unbounded A in a Banach space setting consult [H Chapter II, Section 3].) 

Since the topic of this article is not the behaviour of a single semigruop, 
but the sensitivity of their behaviour with respect to small perturbations of a 
parameter 0, we do not only consider one generator A or one semigroup U , 
but consider parametric families {Ag)g^Q and {Ug)e^Q of generators Aq and 
semigroups Ue = (e''^'')tG[o,oo), with S 9 C R. (Note that 9 is an addi- 
tional parameter and must not be confused with the parameter t of a single one 
parameter semigroup.) 

We further consider a second linear space F that is in duality with E and 
parametric families {TTg)g^Q such that tt^ e F and Agire = for an adjoint 
(dual) a; : F F of Ag. 

We first show in Lemma [1.51 that differentiability oi i-^ Agirg at 9 — is 
equivalent to differentiability of i-)- A^ttq at 9 = provided that 

[A* ~ A*g][TTe - M ^ for 9^0, 
with all limits taken in the weak sense with respect to the given duality. Further 

dA'^TTe I _ dA*gT:Q 



=0 — — |9=o, 



(1) 



d9 ' " d9 

with the involved limits again taken in the weak sense. 

Thus Lemma lTTSl makes the calculation of — -^^^le^o in some cases easier than 
the calculation of ^^\e=o- This fact is exemplified by Remark 14.41 Proposition 
14.51 and Remark 14. 61 in the case that Ag is the generator of the Wright-Fisher 
diffusions. 

Next we prove Theorem 12 .61 the main result of the article. In the case that F 
is a subspace of the algebraic dual E' of E, the theorem provides a formula (see 
Remark [277| for ^^^^^^^no\g=o with [Ug{t)]* denoting the uniquely determined 
adjoint of Ug{t), i.e., a formula for the sensitivity of 1 1-^ [Ue{t)]*7To with respect 
to the parameter 9. This formula involves v := le^o and operators Vo{t) 

given by the series expansion Vo{t) = J2iLi 

It is a simple fact that in the case that F C E' we have [Ug{t)]*TTg — irg, i.e., 
■Kg is a stationary vector of [Ug{t)]* . 

One of the hypotheses of Theorem l2.6l is that E ~ Uje j -^j with Ej Banach 
spaces with respect to the norms such that Ag{Ej) C Ej for all £ 

and the restriction Ag\Ej of Ag to Ej is bounded with respect to It is 

essential for the proof of Theorem 12.61 that the Ej do not depend on 9. Note 
that this hypothesis is fairly restrictive. It prevents us, for example, from the 



^Mappings U from [0, oo) to the space of linear operators U{t) : E ^ E on some linear 
space E such that U{t + s) = U(t)U(s) 
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investigation of diffusion equations with arbitrary coefficient functions. It allows 
however quite interesting insights in the following situation: 

Our abstract results are applicable to the case that E equals the space of 
polynomials on some appropriate real interval and the generators Ag are of the 
form 

Ae:=j2p^{x)-q,{e)-^ (2) 

i—l 

with Pi polynomials of degree less than i and qi differentiable functions. This is 
due to the fact that the operators Ag leave for any k G N the spaces of polyno- 
mials of degree less than k invariant. We apply our results to two examples of 
diffusions operators Ag that fulfill ([2]), i.e. to differential operators of the form 
(HI with n = 2. 

In the first example we demonstrate the applicability of our results to a para- 
metric family of Ornstein-Uhlenbeck semigroups 1 1-^ [Ug{t)]* corresponding to 
generators 

considering without loss of generality the case cr^ = 1. 

OU-semigroups and stochastic processes corresponding to these semigroups 
are frequently used in interest rate modeling. The parameter 9 is interpreted as 
the interest rate to which the process reverts. (Compare with [Ml Vol 2, Chapter 
46] and [HI Section 9.3].) In this example all derivatives can be represented by 
functions. Further it is possible to calculate the derivatives directly since the 
evolution of the OU-semigroup is explicitly given by ([5T|) . Thus the example of 
the OU-semigroup is just of an illustrative nature that does not really rely on 
the developed theory. This is quite different for our second example: 

In our second example we consider a parametric family of Wright-Fisher 
diffusions with mutation and without selection, that can be described by the 
semigroups [Ug{t)]* corresponding to the generators 

^,:=(l-.)^|:+..|+.(l-.)i^. (3) 

Wright-Fisher diffusions are useful tools in population genetics, describing the 
distributions of allele- frequencies in a population (see [5] and Remark l4.7p . Note 
that the stationary distribution ttq of [Uo{t)]* is in the degenerate case 9 = 
and /X > given by the Dirac measure at 0. We calculate the sensitivity ^^\g=o 
of TTe at 6* = as well as the sensitivities Aq^^^ and ^'^g'^" le^o- We obtain 
(Proposition l4.61 formula (|37l) ) that A*^^\g^Q = i.e., the sensitivities under 
consideration are in general not representable by functions or measures on [0, 1], 
but are general linear functionals on the space of polynomials on [0, 1]. 

What makes the concrete calculation of these sensitivities difficult is the 
fact that the involved operators are not diagonalizable. It is however possible 
to construct a basis (see Remark 14. 9|) of the space of polynomials such that ^0 
is almost diagonalizable in the sense that equation (|48l) holds. This enables us 
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to provide a relatively simple recursive formula for the sensitivity -Q^'KQ\e=o 
(Theorem 14. lOp with respect to this basis. 

Although our methods are purely functional analytic and non-probabilistic 
in nature we can — in the case of diffusion processes — interpret the action of the 
sensitivities on the space of polynomials in a probabilistic manner: The action of 
the derivative on the n-th monomial is simply the derivative of the n-th moment 
of the parametric family of probability measures under consideration. This can 
be further interpreted in the case of the Wright-Fisher diffusion (compare with 
Remark [iJ]). 

Derivatives of diffusion semigruops with respect to an additional parameter 

have been dealt with in the context of mathematical finance mainly in the 
context of the stochastic calculus of variations, but also in a PDE context. For 
an introduction to such results consult [9] (especially [9l Theorems 2.2 and 2.3]) 
and [Uj . (For an elementary approach to the relationship of diffusion processes 
and diffusion equations consult [l^)- Derivatives of Markov kernels have been 
considered in |10| and [5] . An extension to derivatives of general operators in a 
Banach space context, relating the derivatives of the operators to the derivatives 
of their stationary vectors has been given in [12] . 

1 Generators and stationary vectors 

Remark 1.1 We follow in the style of presentation of our general functional 
analytic results [71 Sections 16 and 21]. 

Definition 1.2 We say that {E, F, (.|.)) is a dual pairing of the linear spaces E 
and if (.|.) : E x F ^ R is bilinear. We denote by w{E, F) and w{F, E) the 
weak topologies induced by the families of mappings | A* G -F} and 

{/X I ^ S respectively. We say that the dual pairing is separating 

if w{E,F) and w{F,E) are Hausdorff. In the case of a separating dual pairing 
we may identify F with a subspace of the algebraic dual E' and vice versa 
E with a subspace of F' . We denote the spaces of w(i?, F)-continuous linear 
transformations A : E ^ E by Cu,{E) and the space of ■u;(F, £')-continuous 
linear transformation A : F ^ F hj Cw{F), respectively. We say that the 
linear transformations S : E E and T : F F are dual if for arbitrary ^ G E 
and fi G F we have {S^\^i) = (^|T/x). In the case that the pairing is separating, 
we call a dual transformations an adjoint and note that the adjoint is uniquely 
determined. 

Remark 1.3 Let {E, F, {.\.)) be a dual pairing. Let E' denote the algebraic 
dual of E, i.e, the space of all linear functionals (continuous or not) on E. 
Given a parametric family {iJie)eeQ € F^ such that 

e -E {i\pLg - Ho) exists, (4) 
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we let ^&-\e=o G E' denote the unique linear functional such that 



lim e ^(^1/10 - ^io) 

B — ^0 



de 



=oiO- 



(5) 



We call ^^\e=o the £'-derivative oi 9 at ^? = and say that some v £ F 



represents ^j-|e=o if 



de 



-oiO- 



(6) 



In the case that the dual pairing is separating the representative v a F is unique 
(if it exists). 



Proposition 1.4 A mapping S : E ^ E possesses a dual T : F F if and 
only if S £ Cw{E). Further T € Cw{F) (since T possesses the dual S). 

Proof: See |7j 21.1. 

Lemma 1.5 Let {E, F, {.\.)) be a dual pairing. Let G 6 C M with an 
accumulation point of Q. For each 9 d Q let Ag d Cm{E) and irg G F. Denote 
by Ag a dual (the adjoint) of Ag. (Note that the existence of Ag is granted by 
Proposition \1.4\ l- Suppose that: 



and 



(V^gG) A*gTrg=0, 



(V^eE) \ime-\^\[Ag~Ao]*{ng~7ro)}^0. 

6 — >0 



(7) 



(8) 



Then 6 i— )■ AqTTq possesses an E-derivative at that is represented by —u G F, 

%. 6. 

(V^Gi?) limri(^|A5(7r9-7ro)) = (e|-;^), (9) 
g — 

if and only if 9 i-^ AqTTo possesses an E-derivative at represented by v £ F , 
i.e., 

(Ve G E) lim^-i(e|[A; - Al]^^) = 
g — ^0 

Proof: Let ^ e i? be arbitrary. Calculation gives: 
A*g-Kg — A*qt:o 



(10) 



[A*g - Al]Tlo 



I 



(11) 



The limit ^ on the left hand side of equation (|lll) exists and equals by 
^ and ([8]). Thus the same is true for the right hand side, i.e., 



hm ( ( e 



Wg - ^51^0 
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and thus further that 




e 



) 




e 



(12) 



in the sense that if the hmit on one side of equation exists, then so does 
the hmit on the other one. Since G E was arbitrarily chosen ([T2|) cstabhshes 
the equivalence of ^ and (|T0| . 

Remark 1.6 If limg-j-o |(Cl'''6i ^ ttq) = (CKo) for some ttq S F, then © is 
fulfilled with v — AqTt'q. This follows from the continuity of Aq (Proposition 

[El. 

Remark 1.7 If the spaces E and F in Lemma 11.51 are additionally endowed 
with norms and \\-\\f respectively, such that (.|.) is continuous with respect 
to these norms then dH) is implied by 



that is further implied by the norm-Lipschitz continuity oi 9 t-^ Ag and 9 t-^ irg 
at 9 — 0. (Compare with [12] Proof of Theorem 3.1) 

2 Sensitivity analysis of semigroups 

Remark 2.1 Let {H, ||.||) denote a normed linear space and let ||.||-lim„^.oo 
denote the limit with respect to || .|| of the sequence (^n)neN € H^- We say that a 
linear operator A : H ^ H is || .|l-bounded if ||yl||£ := sup^g^^^oj ^j^jp < oo. We 
denote the space of ||.||-bounded operators by £{H) and note that {£{H), \\.\\c) 
formes a normed algebra. If convenient we write ||A|| instead of |1 A||£. 

Lemma 2.2 Suppose that E Uje/-^i linear space and that {Ej,\\.\\j) 
are (for j € J) complete normed spaces. Suppose further that Ag : E ^ E is 
linear, Ag{Ej) C Ej and Ag\Ej G C{Ej). Then for any j G J and any e Ej 



(Ve e E) lim 9-^\\[Ag - Ao]\\ ■ Wng - no\\ = 0, 



(13) 



N 



{tAg} 



71 



:-l|.|lr E 

n=0 



■0 G 



(14) 




ij e E, 



exist, i.e. ^14^ well-defines operators Ug{t),Vg{t) : E ^ E. 
Proof: The Chauchy sequences 




(15) 



converge by completeness of Ej with respect to 
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Remark 2.3 Note that Vo{t)£_ = J*^^ Uo{s)^ds, with the mtegral either taken 
in the sense of Riemann or Lebesgue. 



Notation 2.4 In the situation of Lenima l2.2l we write \\A\\j instead of | 



J- 



Remark 2.5 The following Theorem 12.61 is the main result of this article. Its 
reformulation Corollary 12.71 is concerned with the sensitivity of [Ue{t)]*TTQ with 
respect to the parameter 9 at 9 = 0. Note that [Uo{t)]*TTo = ttq, i.e., ttq is an 
equilibrium for the dynamics governed by t i-^ [Ug{t)]* and thus Corollarv 12.71 
provides formulas for the calculation of the first order effect of small pertur- 
bations of the parameter 9 to the systems dynamics at an equilibrium (of the 
unperturbed system). 

Theorem 2.6 Suppose that the hypotheses o f Lemma \1.5\ and the hypothesis of 
Lemma \2.2\ are fulfilled. Suppose that for any j £ J the mapping 



9^\\Ae-Ao\\j (16) 
is Lipschitz continuous at 9 — and that 

sup |^<c, <oo, (17) 
jG-Ej\{o} II^IIj 

for appropriate constants Cj i.e., ^ i— > (Cl'^'o) defines (for any j £ J) a 
continuous linear functional on Ej. Let Ue{t),Ve{t) : E ^ E denote the opera- 
tors defined by {14^ . Then: 

... The mappings 9 ^ \\Ue[t) - C/o(i)||j and 9 ^ \\Ve{t) - Vb(i)||j 
are — for any j G J — continuous at 9 = 0, 

(ii) y^eE \img^o0-mUeit)-Uo{tm7:o) = -{VomW) 
for any v £ F fulfilling ilO\) . 

Corollary 2.7 In the case that F = E' and '.— /i(C); W6 may reformulate 
conclusion (ii) of Theorem \2.6\ (using Remark ] 1.3\) as 



'^'''^'^^*''''' = lmeit)Uo)=-{VomW) (18) 



89 I 89 
with \Uo{t)]* the adjoint ofUg{t). 

Proof of Theoreml2.6l To prove (i) we just show conitnuity of 6* i-> ||V6i(f) — 
Vb(i)||i at = 0, since continuity of i-> || [U0{t) — C/o(i)||j is proved completely 
analoguous. Let 

(19) 

'^-^ n (n — ly. 
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Lipschitz continuity of 6* |1 Ag— AqUj at = implies that 9 i-^ \\Ag\\j < c < oo 
for some appropriate c > in a 0-ncighborhood (— p, p) Ci &, thus that 



i=0 



<n-c" 



\Ae-Ao\\, (20) 



on (— /9, p) n O and thus further that 



^ \\Ve^N{t)-Vo,N{t)\\j 



N 



71=1 



(21) 



is continuous at 6* = 0. Assertion (i) follows easily from the fact that the limit 
of a uniformly convergent sequence of function that are all continuous at 6 — 
is again continuous at 9 = and that 9 h- !■ ||Vei^Ar(t) — Vo,Ar)(i)||j converges uni- 
formly on (-p,p) to 6* \\Veit) - Vo(t)||j. 



To prove (ii) note that AqTTo — implies 

- {A;r]7To ^ {A;r-'[Al ^ A*]7r,. 

Thus 



(22) 



{[Ueit) - UomM - lim ( J2 ^[(^«)" - (^o)"]eko 

AT— >oo \ n.' 



N 



N 



( ^1 E Sk^^)" - (^5)"]-o ) ihn ( CI E 



(b) 



N 



t [tA, 



*\n— 1 



JV-s-oo 



N 



hm ( [Ae - Ao] V - 



N^oo \ ^ n [n — 1)! 

\ 71=1 ^ ' 



— ^ n (n — 1)! 



(23) 



with (a) and (c) consequences of the || .||-continuity of ^ i-> (CIttq) (Hypotheisis 
(fTT])) and (b) a consequence of (f22|) . 

Let 6'^^[Ag — AqJttq. From (|17p and by the Lipschitz continuity of 

9 ^ \\Ao - AqW at 9 ^ we obtain that 3ij > such that V6l e 6 \ {0} and 

iCWe) - e-\[Ae - Aoieko) < c,||[A, - AoU < <^M\\j- (24) 
From dm) we get G 9 \ {0} and £ Ej that 

{[Ve{t) - VomWe) < iA\\ye{t) - Voimh- (25) 
From ((25)) and (i) we obtain that 
Vj e JVe, e lim{[Ve{t)-VomM < 6j hm || [V^e(i)-K)(t)]e||, = 0. (26) 
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From the weak differentiability of i— > A^tto at 6* = 0, i.e., from ([TU|. we obtain 
that 

V^e-B \im{Vo{t)^We-J^) =0 (27) 
Using (|26p and (P7)) we obtain for ^ ^ E, i.e. for ^ e for appropriate j, that 

lim(Fe(t)ei^9) - (VomW) = lim([l^,(i) - V.mWe) + lim (T/o(t)ei N " = 0. 

(28) 

We finally obtain using and (gS]) that 



hm ^^-l([t/,(^) - UomWo) ^ lim^-i(F,(i)^|[A; - Al]no) 



Thus (ii) has been proved. 



(29) 



9^0 



3 The OU-semigroup 

Definition 3.1 Let / be some interval in M. Denote by TZ{I) the vector space 
of all polynomial functions ^ : / — )■ R. Denote by 71(1)' the space of all lin- 
ear functionals with values in R on Ti{I) and let (.|.) denote the natural dual 
pairing between TZ{I) and ??,(/)'. Note that any such functional is uniquely 
determined on the space of monomials and thus any such functional may be 
uniquely represented as F{^) = { ^"fjf U^o^(^))- 

Example 3.2 Let : 7^(R) 7^(R) be given by 



and let A*g : 7^(R)' 7^(R)' denote the dual of ^e- Let ttq £ 7^(R)' be imphcitly 
given by 

Ve e n(R) (eke) = / e(x)^e("-^)' dx. 



I-K 



Define the Ornstein-Uhlenbeck semigroup 1 1— ^ [f7e(i)]* as the semigroup of the 
adjoints [[/^(t)]* of the operators Ue{€) — e*^^ . Then A^TTg — and the action 
of the Ornstein-Uhlenbeck semigroup on ttq is given by 



/OO 1 
^(a;)^e(--(i-^"')«)' dx. (31) 
-OO 



Further ^e'^ ^'^ and ^e'^ '^"'^ '•'^•'^ are the densities of the normal distri- 
butions N{9, 1/2) and iV((l - e-*)^, 1/2), respectively. Thus by Theorem [53] 
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(ii) and ^ 



mVoitWiy) = (VomW) = lim - Uoimino) 

u — X) 

/OO 1 

e(x)-7r-i/2(e-*-l)a;e^'da;, 



i.e., [Vb(t)]*i^ is represented by the function 7r~^/^(e~* — l)xe^ . It is possible to 
perform the above calculation since we can, in the case of the OU-semigroup, 
calculate [Ue{t)]* and thus [Ug{t)]*TTQ in closed form. Another possibility to 
calculate {Vo{t)(^\i') would be to calculate and to use Remark 12.31 However 
to do this it is again necessary to calculate [Uo{s)]* in closed form. For the 
Wright-Fisher diffusion this has only been achieved in some special cases [5], 

12]. 



4 The Wright-Fisher diffusion 

We intend — in the case of the Wright-Fisher diffusion — to utilize Theorem 
(ii) for the series expansion of limg^oG^^ {[Ueit) — Uo{t)]£^\no) via the series 
expansion of Vo(t). This is done in the next section. 

Remark 4.1 Let k > be fixed throughout this section. For > we define 
operators Ag : 7^([0, 1]) 7^([0, 1]) by 

A.:=(l-.)^|--..A+,(i_,)^. (32) 

Let TTg e TZ{[0, 1])' be implicitly defined by 

e 7^([0, 1]) {(M := J e(x)^|^x''-i(l - xT'' dx (33) 

for 9>0 and by (^Itto) := ^(0). Then A*gng = 0, i.e., 

Vee7^([0,l]) {AeC\7re)=0. (34) 

This is almost trivial in the case that 6 = 0. For 6 > and polynomials ^ of the 
form 

^{x) =x^{l-x)^ -pix) 
— with p{x) an arbitrary polynomial — we obtain (1341) by partial integration 



{Ae^M = 



-^x'-^l - xYO + ^x\l ~ xT-^n + ^x^il - xr 
ox ox ox 



dx = 0, 



which can be further extended to arbitrary polynomials by approximation ar- 
guments. (See also [5] Chapter 4]) 
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Definition 4.2 Let [Ug{t)]* be the adjoint of Us{t) = e*^" with Ag given by 
((32|) . We call the semigroups 1 [Ug{t)]* Wright-Fisher diffusions. 



Remark 4.3 Since for 6 > the function x i-> i^T^ir^a;^ ^(1 — x)'^ ^ defined 

r(e)r(K) ^ / 

on [0, 1] is the density of a Beta distribution, we obtain that {IIttq) — 1 and thus 
further that 9-^{l\TTe - ttq) = 0. 

Remark 4.4 For n > 1 we obtain from ([55)1 that 6 > implies 

1"^^ = r(g)r(^ + . + n) • = n ^T^' (^^^ 



while (a;"|7ro) = 0" = 0. Thus 

\iin0~Ux"\ng-Tro) = lim e-^(x"\Trg) = lim ^ TT 



i— 



Proposition 4.5 Lei &e (7i?;en by i3S\) . Then 



d 
d9 



(C|A;^o)|e=o = lim([Ae - Aol^l^o) = ^(1 - x)^^{x)\no'j - ^'(0). (36) 



Remark 4.6 From ProDOsition l4.5[ Lemma [TTSl and Remark ll.Gl we obtain that 
for Ae given by ([32]) 



d_ 



limJAoC|0-i(^e-7ro)) =C'(0) i.e. := lim^ ^^^-^(Tre - ttq) = — U=o (37) 



and thus further from Remark 12.71 that 



90 '""V ^--v^'^'-' dx' 

Remark 4.7 Calling an element /i e 7?.(M)' a probability-distribution if = 
1 and l/i) > for all ^ > 0, we obtain the following interpretation of our 
Wright-Fisher diffusions t H> [Ue{t)]*: 

Suppose that we start at time in a probability-distribution /i on [0, 1] giving 
us the proportion of individuals — in a large haploid population — that carries an 
allele A. Suppose further that we interpret the parameter k as the mutation 
rate at which allele A transforms into another allele B and 9 as the mutation 
rate at which allele B transforms back into A. Then the probability-distribution 
\Ug{t)\* fjL gives us the proportion of individuals carrying allele A at time t. 
Further the n-th moment {x'^\\Ug{t)]* ji) of \Ug{t)]*ii gives us the probability 
that n individuals independently chosen from the population at time t all carry 
allele A. The probability-distributions ■ng are the equilibrium distributions for 
the respective mutation rates. In the case that 6 = Q and k > none of 
the individuals carries allele A in the equilibrium ttq. Suppose now that we 
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start in the equilibrium ttq, but that the mutation rate 9 is greater than 0. 
Then the probabihty-distribution describing the population evolves according to 
t ^ \U0{i)\*-KQ, and t ^ (a;"|[C/e(t)]*7ro) gives us the evolution of the probability 
that n individuals chosen at random from the population all carry allele A. An 
approximation of the probability (a;"|[J7e(i)]*7ro) for fixed t and small values of 
9 can be obtained by the first order expansion 



(x"|[C/e(i)]*7ro) « (x"|[C/o(i)]*7ro) 



d_ 



(39) 



Since (a;"|[[/o(t)]*7ro) = for > 1 and {x''\[Uo{t)]* ttq) = 1 for n = 0, it sufhces 
by ([55)1 to calculate 



d_ 

dx 



d_ 

09 



[Ue{t)]*TTo 



(40) 



to determine the approximation p9l) . This is done for n = 0, 1, 2 in the following 
example. 

Example 4.8 We calculate the derivative of the 0*'*, I''* and 2"'' moments of 
9 ^ [Ue{t)]*Tro at 6* = by calculating ^[Vb(t)a;*]U=o for i G 0, 1, 2. 



Clearly 



AqI ^ 0,Aox = -Kx and Aqx'^ ^ {-2k ~ 2)x'^ + 2x 



From this we obtain for fc > 1 (by induction on k) that AqI ~ 0, AqX — {—k)^x 
and 



fc-i 



A^x^ = (-2k - 2fx^ + 2 • Y^{-2k - 2)'=-(^+i) {-k^x 

i=0 

and thus further (note that A" ~ id) that 

k=l 



OO 

[Vo{t)]x = Y^ 



fc! 



A 



fe-i 



l = t, 



k=l 



fk 



(41) 
(42) 



k=l 
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and 



oo 

kl 



k=l 



k-1 



i=0 



k-1 



ki 

k=l 
oo ,k 



2k + 2 



(43) 



p(-2.-2)t^.2 ^ 2 ^ ^ M-2At-2)'-(l-(A./(2A., + 2))^) ^ 



g(-2K-2)t^2 



A:! (-2k-2)(1- (k/(2k + 2)) 

2 



K + 2 



1 X. 



From (jH]), (j42|) and (gH]), wc obtain that 



= 0, |:[V,(01=: = and 



-Kt / -(K + 2)t 



1 



(44) 



Remark 4.9 Of course we can also calculate the derivatives of higher moments 
of I— >■ [Ue{t)]*7To at 6* = with increasing computational effort. Moreover there 
exists a basis of the space of polynomials — consisting of the vectors 1, x and the 
vectors ^„ defined in (j46l) below — for that a simple recursion for the calculation 
of \ime^(){9~^[Ue{t) - C/o(i)]Cnko) can be given. 

Theorem 4.10 Let 6„.o,k G M &e arbitrary, let jn.n = 1 arid let for n > 2, 
k > 1 and 2 < m < n — 1 



m(m — 1) 



n{—K, — n + 1) — (m — 1)(— k — m + 2) 

6 



• In 



m=2 



6n,fc := -K6„,fc-i + (n(-K - n + 1)) • 27„,2 



r/ien for n>2 and fc > 1 



^o'K« + K,qx + a] = (n(-K - n + 1))''^„ + 6„_feX 



(45) 

(46) 

(47) 

(48) 
(49) 
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and for bn.o — 



= - \TT^I^\ ^^M,- l_ \ _ X 



lim(e-i[C/eW-C/oW]en|7ro> (50) 

Proof: To prove (|48p we show for = {n{—K — n + 1))'^^^ and i > I that 

Ao[/3^n + bn,i-ix + a\= n{-K -n + + bn,iX. (51) 

This is done by a comparison of coefficients. The coefficient a„ of in Ao[f3£^n + 
bn^k-ix + a] equals n(— k — rt + l)/3as the foUowing calculation shows: 



^I3x" - Kx-^^x'' = n(-K-n+l)px''. (52) 



For 3 < m < n the coefficient a„i-i of ^ in Ao[/36i + box + a] is given by 
n{—K — n + l)/3"fn,m~i as can be seen by the following calculation: 

ml ^^m 2^^ ml ^ml 

Q^m — 1*^ '^T^ n"*^ /^Tn/m n"^ X l3^n.7n~l 

ox'^ ox ox 

= {m{m - l)/37„,m + (m - 1)(-k - m + 2)/37„,„_i)a;™~^ 

= n(-K - n + l)/37„,m_i a;™"\ 

(53) 

with the last equality a consequence of (|45p . Thus it remains to calculate the 
coefficients ai of a;^ and ao of x° = 1. We obtain = since the first and 
second order derivatives applied to the constant function gives 0. We further 
obtain ai = —nbn.i-i + f52"/n.2 ~ bn.i by the following calculation: 

aiX = X-^X^(i"fn,2 - x^-^xl3bn,i-i - KX — x(3bn,i~l j-g^^ 
= [P ■ 27,1,22: + - Kbn^i_i)x = bn^iX, 



with the last equality a consequence of (|47p From ([52|) . (|53|) and (|54p we obtain 
(|5ip and by recursion over (|5ip (with z ranging from 1 to /c) we obtain 
Equation (|49p is a consequence of (|48)) and Remark 14.61 



Finally we obtain ([50| by the following calculation 



u — >-U 

°° J-k 



W , /v^t(iAo)^- 



lim ( - 



N^oo \ ^ A; (fc - 1)! 



/ fe=l 



(55) 

with (a), (b) and (c) consequences of Theorem l2.6l (ii). Lemma [2?2] and equation 
with 6„.o = a = 0, respectively. 
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